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Earlier described model amplitudes are used in this paper to evaluate both 
cross sections and density matrices for two-photon mediated resonance pro- 
^ , duction in e^e~ collisions. All 25 qq low-lying ^5*0, ^Pj and resonances 

\Q \ can thus be treated. Two independent methods are described to obtain the 

■ resonance production density matrices and cross sections. These density 

^ ■ matrices combined with a resonance decay density matrix give the detailed 

angular distributions of the resonance decay products. For two particular 
^ . decays, Xc2^ Xci iJ/'^ the details are given. Several numerical results are 

^ \ presented as well. 



1 Introduction 

When performing electron-positron collision experiments, there will always be data orig- 
inating from (virtual) two-photon collisions. In the latter process single resonance pro- 
duction occurs. When present and future experiments collect more and more statistics it 
becomes worthwile to have theoretical predictions available for observables, which gradu- 
ally become experimentally accessible. 

In a previous paper hereafter called I, a specific class of such observables namely 
form factors were studied. In this paper other observables will be focussed upon, such as 
azimuthal dependences in cross sections and resonance production density matrices. The 
latter quantity combined with a resonance decay density matrix will give the distribution 
of the decay products of the resonance. 

In order to study the above observables one needs a model for resonance production. 
In I such a model was described starting from heavy quarks in a non-relativistic bound 
state. The predicted two-photon decay widths of resonances R turn out to be in reasonable 
agreement with experiment. Once certain modifications are made, also the widths for 
light quark mesons become reasonable. Therefore the model of I is chosen as basis to 
predict observables for the 25 low-lying ^5*0, ^Pj and resonances. In paper I analytic 
amplitudes for 77 — >■ are presented and expressions for form factors are derived. A 
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number of numerical results for cross sections and dependences can also be found 
there. 

In the present paper numerical studies of azimuthal distributions will be given. Their 
characteristic features can be understood from various analytic expressions. Whereas this 
discussion is a simple extension of the results of I, the evaluation of production density 
matrices requires substantial changes in calculational techniques. Moreover, rotations are 
required between reference frames favoured by theory and experiment. A discussion of 
two different evaluation methods will be given in this paper. Finally, in order to make a 
link to the angular distribution of resonance decay products, also one specific decay mode 
for two particular resonances will be discussed. 

Summarizing, the purpose of the present paper is to discuss theoretical evaluation 
methods for resonance production density matrices and to present results thereof. Thus 
the full predictions of the resonance production model of I are in this way completed. 

The paper is organized as follows. Section 2 gives the relevant amplitudes for reso- 
nance production. The direct method of the evaluation of the observables is discussed in 
section 3, whereas section 4 gives an extension of the BGMS |@] formalism required for the 
evaluation of density matrices. After discussing a specific decay model for two resonances 
in section 5, numerical results for several quantities are presented in section 6. 



2 The starting point 

In this section the basic ingredients of the evaluation are listed, such as kinematics and 
matrix elements following from the model of I. 

The two-photon mediated production of a resonance R in electron-positron collisions 
is the reaction 

e+(pi) + e~{p2) ^ e+(K) + e~ip'^) + ^{h) + 7(/c2) 
^e+ip[) + e-ip',)+R{pR). 

The four-momenta pi = {Ei,pi) and p2 = {E2,P2) correspond to the incoming positron 
and electron respectively^, whereas p[ = {E[,p'i) and P2 = {E2,p'2) are those of the 
outgoing positron and electron. The four-momenta of the intermediate photons ki = 
{uji, ki) and /c2 = {uj2, are related to the external four-momenta by 

ki=Pi- p\. (2) 

These two photons subsequently react to form a resonance with four-momentum pr 

Pr = ki + k2= P1+P2- p'l - P>2- (3) 
The virtuality Qi of an intermediate photon is defined by 

Q\^-kl = -{V^-Vif. (4) 

As the virtual photons are space-like, the sign in makes the virtualities Qi real. 

^ We use the standard metric diag((7'"') = (+1, — 1, — 1, — 1). The totally anti-symmetric Levi-Civita 
tensor is defined by £0123 = +1- 
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The invariant mass of the two-photon system is given by 

W = pI= {h + hf. (5) 
The total available energy y/s follows from 

s = {pi+P2f. (6) 

For a typical two photon event W is only a small fraction of \/s. 
Some additional invariants can be defined 

u = {pi-p'2)\ (8) 
U={P2-P',)\ (9) 



and as equivalents of Q\ 



These invariants satisfy 



t = {Pi-Pi)\ (10) 

t'={p2-p'2f- (11) 



s + s' + t + t' + u + u' = W'^ + 8ml, (12) 



where rrig is the electron mass. 

In the lab-frame we chose the following parametrization 



pi = {Eb, 0, 0, -Pb 

P2 = (Ef„0,0,Pb), 



(13) 



p'l = {E[, \pW sin 6^1 cos 01, \p'i \ sin 6*1 sin0i, — \p'i \ cos6'i), 
P2 = {E'^j Ip'sI sin 6^2 cos 02, \p'2\ sin 612 sin 02, \p'2\ cos 6^2), 
PR = {Er, \pr\ sin6'HCos0H, \pr\ sin 6^^? sin 0r, \pr\ cos9r). 

This means that the z axis is along the direction of the incoming electron. The x and 
y axes are chosen such that the xz-plane is the accelerator ring plane and the y axis is 
perpendicular to that plane and is pointing upwards. In these formulae Ej, is the energy 
of the incoming leptons and Pb their momentum. The choice of parameters in (|13D is such 
that we have defined the 2;-components of the positron with an explicit minus sign. 
This parametrization leads to explicit expressions for the invariants 

s = AEl (14) 
Ql = 2 (EbE[ - Pb\P,\ cose, - ml) ^ 2EbE[{l - cosO,) = AEbE^sm^ , (15) 



W = pml + A{Eb - E[){Eb - E'^) - 2E[E'2 - 2\p\\\p'2\ cos ^12 
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(16) 



The approximations are valid in the hmit Ei, > E[» rrie. In formula (|T6|) the angle 9^ 
is the angle between the outgoing leptons in the lab-frame 

cos^^i2 = — COS 6^1 COS 6*2 — sin^^i sin 6*2 cos(0i — 02 + tt). (17) 

The propagators of the intermediate photons will appear as a product of the photon 
virtualities in the denominator of the matrix elements. This is the most dominant depen- 
dence on the photon virtualities. From (^5|) it follows that the photons are most likely 
radiated with low energy under a small angle with respect to the incoming lepton. 

When the incoming lepton scatters over a large enough angle, so that it is detected in 
one of the (forward) detectors, this lepton is referred to as a tagged lepton. From the de- 
tected lepton the four-momentum, and thus the virtuality, of the associated intermediate 
photon can be reconstructed. 

The two methods of calculation both start from the cross section formula 

^ (2vr)MW(A:i + fc2-P) l ^ , dV, dV, 

^J{pi-P2y-mt 4^l'^l (2vr)32El(2vr)32E^^ ' ^ 



with dF the invariant phase space element for the resonance. The amplitude A4 has the 
structure 

M = ^j^f^M^,, (19) 

with 

jT = jT(Ai, a;) = vxApih'vx'M)^ (20) 

= J2'(A2, A'2) = ux'Mh'^^^M, (21) 

or explicitly in the model described (cf. also [0]) in 10 

2 

MCSo; Ai, A2, a;, A'2) = ^6 [h, ^2, ji(Ai, A',), j2(A2, A'2)] , (22) 

Mi'Po; Ai, A2, a;, A'2) = # ([ji(Ai, a;) ■ j2(A2, A'2)fci ■ k2 

-Ji(Ai, a;) • A;2j2(A2, A'2) ■ h] + h ■ h) (23) 
-Ji(Ai,A'i)-j2(A2,A'2)tt'), 



M{'P,; Ai, A2, A'l, A'2, Xr) = ^ {te [£*(Ar), ji(Ai, A'J, j2(A2, A'^), h] 

+t'e [£*(A«), j2(A2, A'2), ji(Ai, A'l), h]) , 

M('P2; Ai, A2, a;, A'2, Xn) = ^ (A^i • hj.^iX,, A;)j2.(A2, A'^) 

+kl^,k2ujliXl, X[) ■ j2(A2, A'2) 
-^i/.j2i.(A2, A'2)ji(Ai, A'l) ■ k2 
-k2,jU^u A;)j2(A2, A'2) ■ fci) e^f^'^iXn), 



(24) 



(25) 



^In these equations we have introduced the shorthand notation e(p, q, r, s) = £°'^''^Pa<lf3'''fSs 
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Mi'D^; Ai, A2, a;, \'„ \n) = ^e*^^'^i\n)h,h, ^^^^ 

e [ki, ^2, ji(Ai, A'l), j2(A2, A2)] , 

where 

Ci = r?o, C2 = %i/iy, 03 = 276^71, Ci = AV?>Wgi, 05 = 8730^2- (27) 
Here we have introduced 

^ 16e^|7^W(0)|^ 

{s + s' + U + U' -%miy+^^W' ^ ' 

with 7^*^*''(0) the (derivates of the) radial part of the wave function in the origin and Cg 
the fractional quark charge. The actual values are given in I. 

Replacing the currents ji, j2 by £1 and £2, the polarization vectors of a virtual photon, 
will give the matrix element of the two-photon reaction, i.e. the second step in reaction (|1]). 
That type of matrix element can be simpliflied by choosing convenient axes to describe 
the photon polarization. It will be the basis of the method of section 4. In the expressions 
terms of the form Si ■ ki or ji ■ ki are omitted, since they vanish in our calculations. 



3 Direct method of calculation 



In this section one way of evaluating the cross section (p^Sf ) and the related density matrix 
will be discussed. The motivation to deal with another method in the next section is that 
we want to obtain the numerical results for observables in two completely independent 
ways. 

The method of this section is a 2 ^ 3 particle cross section calculation without using 
the two-photon intermediate step. For the differential cross section alone it would be 
sufficient to express in terms of the invariants (|6|)-(pl]). In fact, those have been 

obtained and are given elsewhere 0. It makes a very fast numerical evaluation of the 
cross section possible. It is natural to perform this calculation in the laboratory frame as 
specified in equation ([T3|) . 

For the density matrices it is more practical to use a rest system of the resonance. 
Every experimental event can be rotated and boosted back to the resonance rest frame 
(RRF) by successively rotating the momentum in the direction of the 2;-axis and 
boosting it to rest. Thus the required Lorentz transformation is 



LR 



( 7 — 7/?cos0sin0 — 7/3sin {/)sin^ — 7/5 cos^^ 
cos COS Q sin cos Q — sin Q 
— sin cos 

V ~7/5 7 cos sin Q 7 sin sin Q 7 cos Q j 



(29) 



In order to calculate the average density matrix in the RRF belonging to a cross section 
in a certain kinematical region in the laboratory system one proceeds as follows. In the 
laboratory system one generates weighted events. For this we use the event generator 
GaGaRes, which is described elsewhere M. The event is boosted to the RRF with Clr. 
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For the momenta in the RRF the amphtude ^a^.a is calculated, where Aj are the lepton 
helicites and A the resonance helicity. From these amplitudes the normalized density 
matrix 

Ea, A^.a-Aa^.a' /on^ 
Pav = 2 J-* — ' ^'^'^^ 

is evaluated for this event. By performing the weighted sum of these individual density 
matrices and dividing by the total weight one obtains the average density matrix for the 
production of a resonance. Of course, when events with weight 1 are generated the above 
procedure becomes simpler. Sometimes unnormalized density matrices are needed, i.e. 
just the numerator of (|30D. 

It is clear that a method to evaluate amplitudes is needed, preferably a fast tech- 
nique. For this the Weyl-van der Waerden (WvdW) spinor calculus is a good tool, as has 
been shown in the literature Once Dirac spinors, momenta and polarization vectors 
have been translated into WvdW spinors, the matrix elements (p2D-(p6D become spinorial 
quantities. 

There now are two ways to continue. The amplitudes can be expressed in terms of 
spinorial inner products. Evaluating numerically these products and their combinations 
one obtains the amplitudes as a complex number. Another way is to interpret the spino- 
rial matrix element as the trace of a string of 2 x 2 matrices, which can be calculated 
numerically very fast thanks to efficient matrix multiplication in FORTRAN90. The lat- 
ter method has been used in our numerical calculations, whereas the former turned out 
to be useful as independent check. More about the spinorial translation and procedure 
can be found in Here we only mention one particular detail, relevant for comparisons 
and cross checks. 

Usually the spin-1 polarization vectors for a particle with momentum 

k'^ = {k^, \k\ sin6'cos0, \k\ sin6'sin0, |fc| cos^) = (/c°, |/c|sc0, |/c|sS(/,, \k\c), (31) 
are chosen to be 

= ^ (0, =FCC0 + is<p, =FCS0 - ic^, ±s) , 
^l _ ko f \k\ \ 

for helicities ±1 and 0. The spin-2 polarization tensors then usually are taken as 

e^"'{±l) = ^ (£^(±l)e^(0) + e>'{0)e''{±l)) , (33) 

£^''^(0) = (£^( + l)£''(-l) + 2£^(0)£^(0) + £^(-1)5^(1)) . 

Of course this follows a certain convention and other choices may be made as can be seen 
in the literature. In the implementation of the WvdW formalism the conventions of 
have been used. Unfortunately the polarization vectors and tensors in this reference are 
different 

e'^{k) = ^(0, -CC6 ± is6, -css =F ^c<^, s), 

■ ■ (34) 
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e'"'(±2) =e''(±l)e^(±l] 



e'"'{±l) = ±1 (£^(±1)£'^(0) + e^'{0)e''{±l)) , (35) 
When we label the sets of polarization vectors and tensors of equations (|32D and (|33|) with 



the index A and the polarization vectors and corresponding tensors in equations (|34D and 
(p5|) with an index B, the two sets are related by 



''B,0 



4^(0) 



(36) 



These formulae give the relations between density matrices in the two conventions. When 
an event has a density matrix pA, evaluated with the set A, then pB can be obtained by 
multiplying each matrix element with a certain phase factor. In a compact notation these 
factors are summarized as follows 



PXX',B 

with for the spin-1 density matrix elements 

( ° 



PX\',A, 

(j) 2(f) + Tr 

+ 7r 

(J) + tt) 



(37) 



(38) 



\-(20 + 7r) -( 

For the spin-2 density matrix elements we find analogously the set of phase factors 



A0, 



AA' 







20 30 40 
-0 20 30 
-20 -0 20 
-30-20 -0 
-40-30-20-0 J 



(39) 



4 Extended BGMS formalism 

A long time ago in a classic paper p| it was advocated to write the cross section of reaction 
(|I|) as a two-step process. First virtual photons characterized by a density matrix pi are 
created which in a second step produce the resonance R. In this BGMS formalism the 
cross section for resonance production (ITH) takes in first instance the form 



(40) 



QiQ 



where 



Pi 



2 

Q 



I a ^13 



2 \ViPi +Pi Pi 



(41) 
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The index summation in equation (|^ is now replaced by a summation over photon 
hehcities Aj and A'^ with values ±1, 0. Introducing the quantity 



i(27r)^| 6\h + k,-pn)Y.M;,,,^,^M,,,.,^,,dT, (42) 

Aii 



the cross section can be written as 

„2 



, « 1 AiA; A2A' , d^^d^^ 



It should be noted that Mx'^x'^XiX2 is evaluated in the two-photon rest system, the BGMS 
frame, where the direction of the photon originating from the positron is taken as 2;-axis. 
The polarization vectors of both photons now become very simple as can be seen from 
(P^). Also \r is completely fixed by 

Xr = Ai - A2, (44) 

and the summation in (^) consists of only one term. 

The expression for the unnormalized density matrix of an event is obtained from the 
cross section formula when one introduces another combination ppM. 

Saa'= E Pi''''pi''''Mx,'x,'x.x,. (45) 

-^1) -^25 A'^, A2 

A = Ai — A2 
A' = A'^ — A2 

It should be noted that for a specific A only specific Ai, A2 combinations can contribute. 
Using the explicit form of the photon density matrices where only a restricted number 
of elements is independent and using 

,A2,A_R 

where rjji = 1(— 1) for the "normal" ("abnormal") series, one obtains expressions for 
Saa' in terms of Maj'a2'AiA2- every expression also new quantities are introduced, 
generalizations of the quantities aAB and tab of with the help of 

X^^ = {ki ■ - k\kl. (47) 

The results are 

E2+2+ = E2-2- = p++p++M+_+_ = 47x;;p|+p^+a|^, (48) 



p++p°OM+o+o + p?V^+Mo+o+ - 2|pniP2+°|Mo^+o cos(0) 
2y^(p^p°V5 + P?Vrcr,^-4|pr||p2+Vos(0)r#^), ^ ^ 



8 



Soo = 2p++p++M++++ + pfp^'^Moo^^ + 2|p+-||p+-| cos(20)M___ 
-4|pr||p+Vos(0)Moo++ 

-8|pri|pr|cos(0)r#^^ 



(50) 



The various a and r definitions follow from (^Sp-(^D]). In BGMS the following combina- 
tions are used 

cttt = (Ttt + ^TT = , + M+_+_), (51) 



= 4s + Tts = -^={M++oq + M_oo+). (52) 



4a/ X-y^ 



In the formulae is the angle between the two lepton scattering planes in the BGMS 
frame. 

For the off-diagonal elements we do the same. The results are summarized below 

S+_ = e'^^^ (2|pnip^°|e-^^Mo++o - |pMp^°M_o+o - pT\pr\e~^'^ M^+o-) , (53) 

E+o = le'^' (p^+|p^°|e-*^M+++o - Ip^Ip^^Moo+o 

+ |p+-||pf |M__+oe^^ - |prilpMe-''^M__o+ (54) 



+p?0|p+°|e-*'^Mooo+ - \ptWM++, 
For the spin-2 resonances one has additionally 

= le''^^ (p++|p+°|e"^^M+o+- - |p+°|p^+Mo_+_) , (55) 



E2+0 = -e'* (p^+|pMe-2^^M+++_ - |pri|p^°|e~*'^M, 
+ |pMp++Af-+-) , 



00+- 



(56) 



= te''^^ (-|pM|p+V-^^M_o+- + |pnipMe-'^^Mo++_) , (57) 

S2+2- = e^^^^|pM|p+-|e-2^<^M_++_. (58) 

Besides the dependence, for the off-diagonal elements there is also an overall 0i 
dependence, 0i being the azimuthal angle of the incoming positron in the BGMS frame. 
In analogy with the a and the r terms we can introduce 

X+o = I xo+ = :7^==M++o+, (59) 

e+o = ^0+ = :r4^^ooo+, (60) 
2y^X^^ 

C++ = — (oo = _^Moo+-, (61) 

C+0 = -^=M+o+_, Co+ = :^4^Mo++_. (62) 



2^ 






i 








1 




/^77 




i 




1 Xyry 



The ( functions are only non- vanishing for spin-2 resonances. For all above functions aAB 
, . . ., Cab holds: when Cab ~ My^A^AiAa then Cba ~ Ma^a^AsAi- 
The off-diagonal elements can now be written as 



^+0 



'^/x^,{Mpt'\\pt>~'^VRr^s - - pT\pt\e-''^VRC^sT) , (63) 

^^^^y^ (p++|pf |e-^X+o - \ptV^+o + IpMlpflw^^X+o 
-\pf\\pt-\e~''''VBXo^ + pT\pt\-'%^ - \pf\pt^Xo^) , ^ ^ 



= e^^yx,, (p++|p+°|e-^<^C+o - \pt'\pt^VRCo+) , (65) 
S2+0 = -e''^' (pr\pt-\e-''^C++ - Ipfllpf |e-%o + \pr\pt^mC++) , (66) 



= e3*^^yx;;(-|pM|pf |e-^VC+o + |pnipMe-'%+) , (67) 
S2+2- = e^^^^e-2^^yx;;|p+-||p+-|r/^a^^. (68) 

In these formulae rjji is the phase factor arising in equation (^6|). The ^ functions and Coo 
vanish for the ^Pi resonance (Mqoaa' and Maa'oo vanish). 

The other off-diagonal elements can be obtained by using that the density matrix is 
Hermitian 

Saa' = (69) 

and by using that for the polarization vectors in the BGMS formalism the following 
identity holds 

Saa' = i~l)'+''^-y-x- (70) 

The unnormalized density matrix after integration over the phase space of the outgoing 
leptons, denoted by E*"*, is then given by 



The trace should equal the total cross section. When we use the expressions for the 
diagonal elements we indeed obtain the BGMS expression 



da 



(fci-fc2)^-fcffc| 



4:pt^pt'"crTT +'2\pi P2 Ittt cos{2(j)) 



+2pt^pfasT + ^pfpt'-ars + pfpf^ss (72) 



p+VVos(0)rT5 



dVi d p-'a 



The explicit expressions for quantities cr^s, tab, Xab, Cab and (ab follow from the am- 
plitudes MAi,A2,Afl as given in table 1 of I. 

In order to have some analytic results we repeat and extend form factor definitions of 
I and give tables with their forms. Thus we introduce 

aAB = 6iP' - M')8n' ^^^^^^^''^^^'h ABiJ'), (73) 

TAB = S{P' - M')M'^^^^^^^^n^gAB{J''y (74) 
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J' 

0^ 
0^ 



Jtt 



MB 

Jtt 



f: 



TS 



f. 



ss 







K 



( X+uM'^ \ 



2k 



fM^^/ofc 



K 



(Qh< 

\ 2u 







o^m2 Ql 



K 



K 



' X_' 

„2 



2 [2Q'iQl-v{Q^,+Qi,)? 







9tt 






-2/t^ 





0+ 


2.(^) 


/t(#)'QlQ2(X + Z/M2) 


1+ 




«:(^)'QiQ2(^ + g?)(^ + gi) 


2+ 


« (2Q\Ql-v{Q\+Ql^\^ 
12 V / 


(5)' ^ (i(2g?g^ - KQ? + QD) + - - QD) 


2- 


-2k [^J 






Table 1: T/ie form factors required for the cross section fAB, Qab, (A, B = T, S). 



XAB = 5{P' - M^)87r^ ^^^ + ^][^^^^ \ ab{J'). (75) 

Ub = S{P' - M')8n'^^^^^^^^^^mAB{Jn, (76) 

Cab = 6{P' - M')87r'^^^^^^^^^^nAB{Jn- (77) 

The form factors are given in tables |1| and where the quantities k, v and X are defined 
as 

p = h- k2, (79) 

X = X^^ = u^- kjkl (80) 



It should be stressed that for every event the matrix S can now be evaluated. Although 
it is calculated in the two-photon rest frame, which is a frame where the resonance is at 
rest, it is not the RRF of the previous section. For every event one has to relate the axes by 
a rotation which gives rise to a transformation of the density matrix. The transformation 
matrices are given in appendix A. When one likes to compare the average density matrix 
to the one of the previous section, one should also perform the change (pTj) for every event. 
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k+o 


m+o 


1+ 







2+ 






2- 











n++ 






2+ 


«:^^[2g?g2_^(g2 + g2)] 






2- 












Table 2: T/ie additional form factors required for the off-diagonal elements. kAB, f^AB 
and riAB (A, B = 



With the above procedure we extend the existing BGMS cross section Monte Carlo 
program Galuga [0 to a program which can calculate the density matrix S. We then 
evaluate for every event generated by this program the density matrix S which is then 
transformed to the experimentally relevant density matrix of section 3. In this way Galuga 
can evaluate the wanted density matrix. Moreover, we thus have an independent check 
on the GaGaRes evaluation. 

5 Example of a decay model and density matrix 

In general the produced resonance R will decay into some final state X, so that the 
reaction is 

e+e" ^ e^e-R e+e"X (81) 
The total amplitude for this process can be written as 

M = J2aP,M)Ax,V,^. (82) 

Axj^ describes the two-photon production of a resonance with helicity Xr. describes 
the decay of the resonance with helicity \r into the final state X. The factor ^(P, M) rep- 
resents the propagator of the resonance and numerical factors. The total matrix element 
still depends on the external four-momenta. No implicit integrations over the outgoing 
momenta have been carried out at this stage. The square of the total matrix element is 
given by 

Y.\M\' = UP,M)\' A,A',^^A«Ai, = Up,m)\'i:t{av*). (83) 

The summation on the left hand side represents a summation over the helicities of the 
initial and final state particles. The quantities ^a^a' and 'Dxp.x' are Axp, A\i and Vxj^Vy , 
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summed over the helicities of all particles but the resonance. These are the density 
matrices for the production and decay of the resonance. For a spin-J resonance the 
density matrices are formed by (2J + 1) x (2J + l)-matrices. For actual calculations a 
choice of polarization vectors and tensors has to be made. The results of this paper are 
obtained in the conventions (0) and (P^l). 

For R we first take the spin-2 state Xc2 and consider the specific decay mode 

Xc2{pr) -> l{h)J/^{h). (84) 

Later on also the possibility of the J/ ip decaying into a lepton pair will be discussed. The 
choice of this example is prompted by its experimental relevance. Since we expect that a 
similar Xci decay may cause a contamination of the events (|8^) in an actual experiment 
we shall also discuss the xa decay at the end of this section. 

For the amplitude M. of the decay of the spin-2 resonance we use the amplitude that 
is given in equation ( P3| ) 



M{\ Ai, As) = Q [(fci ■ k2)eT{\i)er{\2) + (£*i(Ai) ■ e^iX^W^k^^ 

-{k, ■ el{X2))er{Xi)k!^ - {k^ ■ e\{X,))er {X2)kl] (^5) 
= C4Fr^(Ai)F2V(A2)e;:(A), 

where F"^ is the field strength 

Ff{X,) = kte^,{X,)-kU'^{X,). (86) 

In these equations the index 1 refers to the photon and 2 refers to the massive J/ip. The 
A's refer to the helicity of the associated particles. The parameter C4 is not specified here, 
but its complex square will contain the width for the decay of the spin-2 particle into 
these two spin-1 particles. It will be related to C4. 

The density matrix T> for the decay of the resonance is then constructed by 

Vxx' = E -^(^' A2)A^*(A', Ai, A2). (87) 

In the evaluation of this density matrix we make use of the invariant spin summations for 
spin-1 particles. For the massless photon this relation reads 

^ei,a(Ai)eJ.^(Ai) = -c/„/3, (88) 

Ai 

whereas for the massive J / ip it reads 

E e,A^,)elpi^2) = ~9a, + (89) 
A2 

where M is its mass. In fact, due to the gauge invariance of the field strength in equation 
(P5|) only the 9^/3 term in equation ( |5^ contributes. This leads to the following expression 
for the density matrix for the decay of the resonance 

-Dxx' = |54p [(fci ■ fc2)£-^(A)4^(A') + 2{k^ ■ £(A) ■ k^){k^ ■ e*{X') ■ k,) 

(90) 

{kl + A{ki-k2)){k,-e{X)-e*{X:)-k,)]. 
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In this equation a polarization tensor in between two dots indicates that one Lorentz 
index has to be contracted with the four-vector on the left side of the tensor whereas the 
other Lorentz index has to be contracted with the four-vector on the right side of the 
tensor, i.e. 

p-A-q = A^^p^q,. (91) 

In the derivation we have used that the final state photon is massless and we have used 
the properties of the polarization tensor to replace ^2 in contractions with it by —ki. 

The density matrix for the decay of the Xc2 is constructed in the rest frame of the 
resonance with the z axis as polarization axis and some x and y axes. The density matrix 
will depend on the polar and azimuthal angles of the photon in this reference system. 
The specific choice of the axes in the Xc2 rest system is at this point arbitrary. In the 
applications they should be the same as for the Xc2 production density matrix. The 
polarization tensor looks very simple and one can use equations (0) and (|35|) to verify 
that the following useful relations hold 



In this frame the four-momentum of particle 1 can be parametrized by 
fcj' = |fc|(l, sin 9 cos (p, sin 9 sin 0, cos^^) = |A;|(1, s cos0, s sin0, c). 



(92) 



(93) 



It can be shown that with these definitions the density matrix for the decay of the 
spin-2 particle can be written in a compact form as the sum of three matrices 



2 

5aa' + Wvxvl, + (1 + «)i?AA' 



In this expression 5 is the Kronecker symbol. The vector v is given by 

^2 1 „2 



V = (t;2,t'i,fo,t'-i,t'-2) 



-e^'t -see''', -^{Se^ - 1), see-'t -e-^'^). 



The tensor R is given by (the indices run from 2 to —2) 



R 



2 


cs i(f> 
2 ^ 


_i_p2i<j) 

2^6^ 


cs —i<j) 
2 ^ 




sc i(l> 
2V^ 


4 


„-2iij) 
2^6^ 


sc „—i4> 
2V6 


3c2+l 
6 





s^ -2i(j) 
4 


sc ~i 
2V6^ 








s2 -2i4 
2v^^ 







s r.2i4 








2V6 2^/W^ 



_cs 14 
2 ^ 



(94) 



(95) 



(96) 



V u u - ~e ~ J 

In the second line of equation (p^) the expansion parameter a has been introduced. 



a 



Ml 



(97) 
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Others [§] use a simplified expression for tlie matrix element, based on an electric 
dipole transition 

-M = Q [{h ■ h)eT{Xi)er{X2) - {h ■ £KAi))A:^r(A2)] e,,{X). (98) 
This leads to the following density matrix 

Vxy = IqP [(^1 ■ hfSxx' - j^ikl + 2{h ■ k2))\k\^vxvl, + {k^ + 2{k^ ■ k2) 

(99) 

2 



In the limit of vanishing a the density matrices of equations (0) and (0) become equal. 
For the actual a, based on M = 3.09687 GeV, Mr = 3.55618 GeV and therefore 
a ~ 0.242, there are some differences, e.g. for equation (|93) 



V22 + V_2-2 ~ 1 + 1.570 ■ + 0.0190 ■ c^ (100) 

and for equation ( ^9]) 

©22 + ^^-2-2 ~ 1 + 1.028 ■ - 0.00956 ■ c^ (101) 
whereas for a = one would have 

V22 + V^2^2 ~ 1 + (102) 
These distributions can be compared to the experimentally observed decay distribution 



Tl| of the Xc2 produced in pp collisions, 

V22 + ^^-2-2 ~ 1 + 1.96 • c2 + 0.0142 . c^ (103) 

where the central values of the coefficients are taken and the errors are omitted. In section 
6 some numerical results for the distributions will be presented. 

One could also take into account the possible decay of the massive J / ip into two leptons 

Xc2{k) ^ j{k,)J/ij{k2) J{k,)l+ip,)r{p2). (104) 



The matrix element for this reaction can be obtained from equation ( ^5]) by replacing the 
polarization vector of the J/ip by the lepton current 

J:^s = MP2)Yvs{Pi), (105) 

where r and s denote the spins of the final state leptons. In the following the lepton mass 
mi will be neglected, which is a reasonable approximation for this decay. 

In the construction of the density matrix for this decay one has to use that for the 
lepton current one finds 

E Jr'sJ^rs" = 4KP2 + " 1^19^ = 2[k^k^2 ' I'l" " ^19^- (106) 

r.s 
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In this equation the vector I = Pi — P2 has been introduced. This has been done as the 
term proportional to k2k2 will not contribute to the density matrix. For / we use the 
following parametrization 

r = {h,l,Jy,Q (107) 

in the above chosen reference frame, where (|93D holds. The density matrix now reads 



+{2kl + l^){ki ■ s{\) ■ ki){ki ■ e*iX') ■ ki) 

-{k, ■ k2)[{k, ■ e{X) ■ k,){l ■ e*{\') + e{X) ■ l){k, ■ e*{X') ■ k,)] 

+ (2(A;i ■ k2) + kl){kr ■ e{\) ■ l){k^ ■ e*{\') ■ I) (108) 

+ [{klf + 4(fci ■ k2)kl + {h ■ lf]{ki ■ e{X) ■ 6*{X') ■ ki) 

+ih ■ k2){ki ■ l)[{ki ■ e{\) ■ e*{\') ■ I) + (l ■ ^(A) ■ £*(A') • ki)] 

+{k,-k2f{l-e{X)-e*{X')-l)). 

In the derivation we have used that (^2 ■ /) = This density matrix can also be written 
in a more compact form 

Vxy = |£p ((fci ■ k2fkl6xy 

+ [2kl + fm^vxvl, - ■ k2)\k\''[vxxl, + x^vl] 
+ {2{ki-k2) + kl)\k\''wxwl, (109) 
+ [{klY + 4(A;i ■ k2)kl + (fci ■ Ifm-'Rxy 
+ {h ■ k2){k, ■ l)\k\[Sxx' + S*y^] + {k, ■ k2)^Txy 
In this density matrix the coefficients can again be written as functions of a 

+ j^wxwl, + (1 + « + cos2 e*)Rxx' (110) 

2A/Vl-a COs6'*(S'aA' + Sy^) + -^Txx') , 

where a convenient expression 

Ml - 

ki-l= cos 9* (111) 

has been used. The angle 6* is the polar angle of the outgoing in the rest system of 
the J/ ip where the z axis is given by the direction of the boost from the RRF to the rest 
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frame of the J/^/^, i.e. opposite to the photon direction. In the expressions some additional 
vectors and tensors have been introduced. 



.sL 



w 



Txx' 





^ + /+c],--^[s(/. 


COS 4> + ly sin 


0) - 2dJ 


1 


+ /- 




[^x + ^j/ ~ 




z 






^ _ s'- pi'l' sh picp 
2 2 


2V6 












Cl- 1 7 

2 2 «^ 


ls + 

2^6 


4 ^ 


n 




Sw — 


si- -iff) Icz 
2a/6 2a/6 


Q L .^y 


] 2a/6 


2a/6 






^e-^-^ 


/* 

2\/6 


2 


d+ 
~2~ 






^ 


s^- ~i(t> 
2V6 


slz Q—id 

2 


2 / 






2V6 








\ 


l-lz 
2 




2V6 


4 









l-lz _ 

2V6 




2^6 


'+ 
2V6 







4 


2V^ 


_l[|/|2 + ;2] 




V 





275 


l-lz 
2 







sl_ 



In the expressions the variables /+ and /_ have been introduced 



(112) 
(113) 



(114) 



(115) 



(116) 



I— Ix ily 

In addition, we have also introduced 

Ixy — ix cos 4> + ly sin 4>, (117) 

l^, = 2d_ - she-'''', (118) 
= 25/^6^"^ - cl+, (119) 

lsz^^e"'' + ch. (120) 

The decay density matrix is now completely specified once the lepton momenta in the 
Xc2 rest system are inserted in 1^. It is sometimes convenient to parametrize 1^ in terms 
of decay angles 9* and 4>* , 



s* — sin 9*, c* — cos 9*, d^ — cos 0*, 
of the in the rest system of the J/ip. 
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sin (f)*, 



(121) 



The components of the four-momentum I in the RRF then become {c^ = cos (p, 



sin0, the angles of equation (|93|)) 



7 _ Mr-M^ 



U = M{s*{c*cc^ + s*s^) - c*sc^), 
ly = M{s*{c*cs^ + s*c^) - c*ss^), 
h = -M{s*cp + c*c). 



(122) 



Since the Xc2 will be predominantly produced in helicity ±2 states, we have a closer 
look at those density matrix elements. In the limit of vanishing a the sum of density 
matrix elements P_2-2 + ^^22 reads 



■2-2 



22 



|g|2 



2 j M^{2 + ^{W^2WU+W2W;) 
+ ]v72(^-2-2 + ^22) + R-2-2 + -R22 J 



(123) 



Inserting the expressions for the vectors and the tensors gives 
^ ^ , ,2 fMl-M^Y 2. 



2M2 



(124) 



Integrating over the angles associated to the outgoing leptons yields for the sum 
T'_2-2 + "^22 a result that is proportional to 1 + c^, which is in agreement with the 



previously found results for the decay of a Xc2 into J /ifj and 7, equation ( |102[ ) in the limit 
a = 0. Using the expressions for and Z_ results in 



2-2 



+ n 



22 



M'^il + c^) ( 1 



1 r 



,s c^c-c s) + is 



■ (125) 



For helicity-2 predictions these results can be compared to the weight function for the 
angular distribution, /|a|=2, presented in and |TT[ 



f\x\=2{e. e\ 0*) = \AI{1 + c*2)(l + 6c2 + c') + AI{1- c*2)(l - 



+ 1^2(1 + c*2)(l - 2c2 + c^) + ^A2Aicl2s*c*s{^ + 3c 
+^yl2Ao(rt2 - sf)s\l - c') - ^A,A^cl2s*c*s^c. 



(126) 



A pure electric dipole transition would result in amplitudes Aq = 0.316, Ai = 0.548 and 
A2 = 0.775. Experimentally lll| the values Aq = 0.21, Ai = 0.49 and A2 = 0.85 have 
been found. It turns out that the compact expression ( |125| ) amounts to expression (|126|) 
with the pure electric dipole values for A^. 

After the above detailed discussion of a Xc2 decay, we now briefly mention a similar 
Xci decay, since it may experimentally contaminate the Xc2 events. Therefore we need a 
model for the decay 

XciiPR)^lih)J/i;ik2), (127) 
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from which one can derive the density matrix. Again, we use the 7*7* amphtude where 
one 7* is replaced by the J/ip. We thus have from (^) 



M{\ Ai, A2) = hM\[e{X), £*(Ai), £*(A2), h]. 
The density matrix now becomes 

iki ■ k2 



1-2- 



M2 



M2 



(e(A)..*(A')) 



When ki is again parametrized as (0), one obtains 

{ki ■ k2 



1-2- 



M2 



ri2 * , (^1 ■ ^2) 

k\ u\Ux H — - 



where 



u 



[mi,Uo,M-iJ 



^e'"^ -c —e 
V2 ' 'v^ 



■i(j> 



:i3r 



Since in the production the helicity ±1 states are favoured, the angular decay distri- 
bution will now be approximately of the form 



1-1 



1 - 



2M2 - Ml 



2M2 Ml 



^c^ ~ 1 -0.216■c^ 



(132) 



which has a behaviour opposite to the decay distribution for the Xc2^ e.g. ( |102| ). When 
forward or backward parts of the angular distribution "^J/ip are experimentally not ac- 
cessible this will affect the Xc2 decay more than the Xci decay. Although the GaGaRes 
predicted two-photon mediated production of Xci ^ is only 15% of that of Xc2 

0], the Xci contribution to the actually seen '^J/ip events may be higher than this 15% 
because of the opposite decay distribution. In the experimental determination of the 
7*7* production through its Xc2 iJ/i^ decay it would be worthwile to estimate the Xci 
contamination of the events. 

Such a contamination would somewhat reduce the really present 7*7* — > Xc2 produc- 
tion and therefore the derived Xc2 77 width. This remark is relevant since at present 
the measured two-photon width is a factor 3 to 4 higher than the determination from pp 
collider experiments . More recently it has been argued that this factor is about 2 to 3 



It is clear from these remarks that a precise measurement of the 7J/V' angular distri- 
bution could clarify the above situation. Tools to calculate various density matrices will 
then be indispensible. 

One may wonder whether the above amplitudes predict the ratios of the widths Fi 
and r2 for the decay processes (|127|) and (|84D correctly. This is considered in section | 



6 Numerical Results 

In this section numerical results for several quantities are given. They have been obtained 
by using the event generator GaGaRes while extended Galuga results served as a 
check. 
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6.1 Azimuthal distributions 



Before focussing on azimuthal distributions it is useful to show the shape of the cross 
section da/dQ with its marked decrease with Q. This is displayed in figure |l|, where 
every curve has an arbitrary normalization. 

Next, results for the A0 = 0i— 02 distributions for the ^5*0 and ^Pj bottomonium states 
are given for y/s = 190 GeV. From figures ^-0d it is seen that there are marked differences 
between distributions for different resonances. The characteristics can be qualitatively 
understood when the dominant features of the cross sections are considered. Similar 
qualitative arguments, but without numerical results have been given in the literature 
||T3| in connection with Pomeron production of resonances. In the first place it is clear 
from figure |I| that the largest contributions to the cross sections come from the region 
Ql < M2. In those regions A0 (cf 0) such that a qualitative understanding of the 
distribution in the BGMS cross section formula ( [72D is sufficient to explain the features of 
figures 0a-^. The behaviour of the BGMS formula for the various resonances now follows 
from the form factors /ab and qab in table |1| on page |Tl] and the small approximation. 
One then obtains distributions, where different photon density matrices play a role. 
Using the approximate numerical relations valid at y/s = 190 GeV and in the small 
region 

p++^p+-, p+°^1.5p++, p°°^2.3p++, (133) 
one arrives at the following shapes for the distributions: 



^ ~ (pVpV - P\-Pr cos(20) 



~ '^^P^P^ sin^ 



(134) 



I ~ 4«: (^) (p^p++ + pt-pf cos(2 
~8«:(^)prp++cos2 0. 



(135) 



g ~ 4«: [^)" (p^+p++ - pt-pt cos(20); 

+4^^ (f )' [pX^p'^QI + pfp^Q? - M'pfQ^Q. cos 0) (136) 

||p++p++ (m2(Q2 + Q2) _ 2M2QiQ2 COS + (Q? - Qi)2 sin^ 



Xh2 ■ 

It is clear that the sin^ and cos^ and constant distributions show up in figures ^ a,b 
and d. A more complex structure arises in figure |^ c, but also here the cos0 and sin^0 
distributions can be recognized on a constant background. 

The bb state has not been included in the plot as it would lie exactly on top of 
the plot for the ^Sq bb state. 

For completeness we note that for the cc states a similar A0 behaviour has been found. 
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6.2 Density matrices for resonance production 

For the cc resonances Xc2 and Xci we have studied how the different hehcities of the 
resonance contribute to the cross section as a function of the Qf cuts. The results for 
the diagonal elements in the RRF are given in figure |^ In these plots both Q\ and Q\ 
have to be greater than the value given on the horizontal axis. From the right plot in 
figure one can see that in the absence of cuts the dominant contribution comes from 
the helicity-2 components. This is in agreement with |jT^ where a similar dominance 
was found for the helicity-2 component for the /2(1525) resonance. For high cuts the 
helicity-0 component starts to dominate. In the BGMS-formalism this contribution 
comes from the ass component that is proportional to Q\Q\. This result also agrees 



with theoretical predictions found by Close |]T5|, in which he states that the helicity-0 
contribution is of the order 0{tt' M"^). From the left plot in figure ^ one can see that in 
the case of the Xci resonance in the absence of cuts the helicity-1 component is dominant. 
It would be interesting to verify this statement with experimental data. For higher Qf 
cuts again the helicity-0 component starts to dominate. 

6.3 Density matrix for the decay Xc2 iJ/'^ 

In section 5 we have constructed the density matrix for the complete decay of a Xc2 
resonance into a, J /ip and a photon, where the J / ip subsequently decays into electrons or 
muons. We can now use equation ( |5B| ) to combine the density matrix for the two-photon 
production of a Xc2 resonance with the density matrix for the complete decay of the 
resonance to get the full matrix element squared. The average production density matrix 
is given in table 3. Again it is clear that the helicity-2 diagonal elements are dominant. 
The combined production and decay matrix elements (|9^) have been used to generate 
the distribution of the angle between the outgoing photon and the boost direction in the 
RRF. The results are given as data points in figure ^j. 

In this figure also results are included which are obtained from the assumption of 



pure helicity-2 production. Under this assumption the V22 + 1^-2-2 distributions ( |100| ), 
( |101|) or (|103|) can be chosen. Comparing the data points of the full calculation and 
distribution ( [LOOj ) shows that the production is indeed helicity-2 dominated. Comparing 
the distributions (|101|) and (|103|) with the model distribution (|100|) shows that the model 
is closer to experiment than the simple dipole distribution. 

6.4 The ratio r(xci ^ lJ/i^)/r{xc2 ^ iJ/i^) 

As the trace of the density matrices for the decay of a Xcj resonance into a photon and 
a. J/ip yields the matrix element squared for this process, we can use equations (p^) and 
( |130|) to calculate the ratio of the two decay widths. We assume C3/C4 = C3/C4, i.e. we 
assume that the replacement of a 7* by a J/-?/; introduces in both matrix elements (p^) 
and (1^) the same scale factor. The mass of the Xcj particle is denoted by Mrj. Using 
the expressions for the traces of ( p^ and ( |130| ) yields for the ratio 

_ r(xci ^7^/^) M\Ml^+M') ( MrX (o^\ 1 ..... 

r(xe2-7^/^) \MnJ \(3) IQ-ha + a^' ^'"'^ 
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A' = +2 


A' = +1 


A' = 


A' = -1 


A' = -2 


A = +2 


0.47 


2.0 ■ 10-^ 
4.72 


0.012 
6.27 


6.0 ■ 10-^ 
5.56 


0.029 
0.13 


A = +l 




0.022 


5.4 • 10"^ 
1.53 


4.1 ■ 10-3 
0.023 


6.0-10-^ 
2.41 


A = 






0.012 


5.4 • 10-^ 
4.67 


0.012 
6.27 


A = -1 








0.022 


2.0 ■ 10"^ 
1.58 


A = -2 










0.47 



Table 3: The normalized density matrix pxy = |pAA'|e"'^^^' for Xc2 production at y/s = 
190 GeV in the absence of cuts on the external particles. In each entry the absolute value 
Ipaa'I of the density matrix (upper value) and the phase (pw (lower value) are given. The 
open boxes are filled by Hermitian conjugation of the other elements. 



where /? has been defined hke a 



2 



M' 

= 1 - ^ ^ 0.225. (139) 
Inserting the numerical values gives 

Rth ^ 0.8969. (140) 

This number should be compared to the experimentally measured |^ ratio 

p rtot(Xci)Br(xci ^iJ/'ip) naa^mr: ^^A^^ 

Rexp=T^, rg-T —— = 0.89 ±0.15, 141 

rtoi(Xc2)Br(xc2 iJ/^) 

where in the calculation of the error the different contributions to this error are taken to 



be independent. Using the data selection of [0] this number becomes 

Rexp = 0.76 ± 0.26. (142) 
These experimental numbers are compatible with the above theoretical prediction. 
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A Transforming the density matrices 

In section 4 we have constructed the complete density matrix in the BGMS frame. How- 
ever from an experimental point of view the RRF is a more convenient reference frame. 
As the elements of the density matrix depend on the chosen polarization vectors/tensors, 
the density matrix is frame dependent. When one wants to rotate the quantization axis 
over the angles Op and 0p the density matrix p changes into the density matrix p according 
to 

Pw = Ax^,Ay^,pf,f,,. (143) 

Thus a density matrix p calculated in the BGMS frame can be transformed to the RRF. 
In the BGMS frame the boost direction which was needed to get to the Xc2 rest system 
is characterized by 9p and (pp. 

For spin-1 resonances this transformation matrix A can be shown to be 



A 



A/1 



^ 1+c —i<t> s l—c id> ^ 
2 ^ 2 ^ 

__S_ 1±C 

V2 2 ^ / 



~73' 

l-C 

2 ^ 



(144) 



The order of the indices X, p is taken to be (+, 0, — ). For spin-2 resonances this trans- 
formation matrix reads (with order +2, +1, 0, —1, —2) 



A 



-2id 



(l+c)s 



2 V 2 



i^-c)s ( 1-c 



-,2i(t> 



(l+c)s 2i<t> {l+c)(2c-l) U /l" (l-c)(2c+l) {l-c)s 2ii 

2 2 y 2 2 2 



1 /3„2„-2j(/) 

2 V 2 



sce~ 



{l-c)s 2i^ {l-c){2c+l) 
2 2 



\ 



1-c 
2 



-2i<j) 



2 

-CS 

£. /3 
2 V 2 



1 3„2p2i4 

2 V 2 



3 (2c- 1) it/, (l+c)s 2i(/. 

2 2 2 



(l±£)£pi</> 
2 



^2i(j} 



(145) 



In the expressions for the transformation matrices we have introduced the shorthand 
notation 



sm 



J p. 



(146) 
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Figure 1: The Q spectrum of the photon radiated by the incoming positron for the generated 
bb states, the rji, (continuous line), the Xbo (dashed line), the Xbi (dotted line) and the Xb2 
(dash-dotted line). The events were generated at y/s — 190 GeV without cuts on the 
eocternal particles. The vertical scale is logarithmic. 
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Figure 2: T/ie A(f) spectra for the bb states rjb (a), Xbo (b), Xbi (c) and Xb2 (d). The events 
have been generated at y/s = 190 GeV in the absence of additional cuts on the external 
particles. The vertical scale is linear. 
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Figure 3: Contributions of the different diagonal density matrix elements of the xa 
resonance (left plot) and the Xc2 resonance (right plot) to the total cross section at 
= 91.5 GeV, for different Qf cuts. (I.e. Ql > 0, 1, 5, . . . , 25 GeV^ .) 
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Figure 4: The distribution of the polar angle of the outgoing photon in the RRF. The data 
points show the result from the calculation using the combined matrix element both for 
the production and decay. The dashed line shows the distribution ( \10(\ ) for the complete 
density matrix for the decay of the resonance where the production is purely helicity-2. 
The dotted line represents the distribution ^101\ ) according to the pure dipole transition 
and the dashed-dotted line shows the distribution that has been experimentally observed 
/|77|/. The vertical scale is linear. 
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